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Abstract 

A generalisation of the non-perturbatively stable solutions of string equations which 
respect the KdV flows, obtained recently for the (2m — 1,2) conformal minimal 
models coupled to two-dimensional quantum gravity, is presented for the (p, q) 
models. These string equations are the most general string equations compatible 
with the q-th generalised KdV flows. They exhibit a close relationship with the bi- 
hamiltonian structure in these hierarchies. The Ising model is studied as a particular 
example, for which a real non-singular numerical solution to the string susceptibility 
is presented. 
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1. Introduction and Conclusions. 

Pure two dimensional quantum gravity is known to have a KdV flow symmetry 
to all orders in genus perturbation theory. In a series of papers [l]-[4] it was 
established that a complete formulation of non-perturbative pure two dimensional 
quantum gravity can be developed from the single principle that this symmetry 
is respected non-perturbatively i.e. that the KdV flows are exact. In particular 
the principle leads to a unique string equation which has as special cases the non- 
perturbatively sick solutions of hermitian matrix models, but also has a uniqueAl 
real non-singular solution. The above discoveries are briefly reviewed in section 2, 
together with the matrix model reasons for expecting these successes. The primary 
reason for the present paper is to show that these successes generalise to 2D quantum 
gravity coupled to a general conformal minimal model. For (p, q) matter the single 
principle is that the q-th generalised KdV flows (known to exist perturbatively[7]) 
are preserved non-perturbatively. We will demonstrate that this leads to unique 
string equations for each system as conjectured in ref.[2] (and indeed for massive 
theories interpolating between all (p, q) critical points for given q) and display a 
uniqueA2 real singularity-free solution for the Ising model coupled to gravity. 

The success of the present formulation demands an analysis of its principle: 
non-perturbative preservation of flows. Let us first note that if this principle is 
discarded then other formulations are possible[ll]; Clearly some input is needed 
to define 2D gravity (a.k.a. ID string theory) beyond the genus expansion. One 
might hope that unitarity ensures a unique non-perturbative extension. For 2D 
string theory (the lowest dimension where the concept of an S'-matrix makes sense) 

Al We believe. In ref.[2] it was proven that there is at most a discrete number of 
such solutions with real asymptotics, and a numerical study uncovered only one. 
A2 Subject to similar caveats. 
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this appears to be insufficient [5]. In this case it is a natural conjecture that the 
perturbation expansion has a form of KP flow symmetry and it would certainly 
be interesting to trace this out and determine whether or not an exact KP flow 
symmetry picks out a unique unitary non-perturbative extension. Of course some 
other general constraint (causality?) conceivably might provide the missing infor- 
mation but this is not really the point: If one succeeds in providing an 'unprincipled' 
extension (in the sense that the perturbation theory must be separately determined) 
then non-perturbative string theory remains logically incomplete. 

Having argued for principles, what form should they take? A common attempt 
in the past has been to try to formulate non-perturbative string theory from some 
symmetry principle on the world-sheet i.e. a symmetry of 2D quantum gravity 
considered on a single (possibly pinched) genus. In view of the work of the last 
two years on low dimensional string theory, and of the simplicity in this case of 
the world-sheet theory when appropriately formulated, this now surely seems as 
unlikely as expecting, say, non-perturbative QCD to arise from some symmetry of 
the world-line. Indeed in a second quantized theory one expects the symmetries to 
be best manifested on the second quantized fields. Recalling that their expectation 
values - the background fields - are nothing but the couplings t r in string theory, 
we expect the symmetry to be manifested as active transformations on the t r . 
Thus a non-perturbative string theory symmetry principle should be a symmetry 
of "theory space" , the space of all world-sheet couplings, and not of a single world- 
sheet theory. This is precisely what the KP (generalized KdV) flows are. Thus 
it seems highly probable that this symmetry principle is a hint of a much larger 
symmetry determining the non-perturbative form of 2D quantum gravity coupled 
to general confer mal matter. 

In this paper our primary purpose, as already mentioned, is to develop the 
formulation for gravity coupled to (p, q) matter. The first steps towards this were 
already taken in ref. [2] where, motivated by a matrix model whose eigen- value space 
was IR-i-, use was made of [P, Q] = Q: the appropriate generalisation of Douglas' 
[P, Q] = 1 formalism[7]. By using the Lax pair formalism for (generalized) KdV 
flows we will see that in this case the Douglas' formalism is nothing but a trivial 
expression of scaling: i.e. it follows directly from the fact that KP flows have a 
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grading. Many other intimate connections with the reductions of the KP hierarchy 
are uncovered. For example the classical VFA((/)-algebra of the second hamiltonian 
structure is seen to play a central role. The integrability of the hierarchy implies 
the existence of a first integral of the scaling equation: our string equation. The 
existence of a natural 'gauge' parameter a - the boundary of 1R+ and physically 
a world-sheet boundary coupling - is seen to be a consequence of coordination 
of the bi-hamiltonian structure. Indeed another reason for the present paper is to 
provide a complete discussion of this parameter, partial results having been reported 
earlier [3] [11] [4]. 

In section 5 we recall[2] that the L_i symmetry of KP flows is not a symmetry 
of the vacuum, leading to an analogy with spontaneous symmetry breaking in which 
a is identified with the Goldstone boson. Since, as we will show, all the WX(k) n 
generators with negative index n are 'spontaneously broken' there are further 
generalizations of the a parameter for the (p, q) models with q > 3. For pure gravity 
a is the boundary cosmological constant. For the Ising model it is the boundary 
magnetic field. There are two other parameters for the Ising model associated with 
W-i and W-2- The latter we tentatively associate with the boundary cosmological 
constant - a parameter missing in previous formulations [19]. 

The structure of the paper is as follows: Section 2 is a short review of previous 
work followed by a review of our construction for the (2m — 1,2) models. We 
emphasise the Lax pair formulation and bring together the earlier results on a. 

Section 3 deals with the Ising model coupled to 2D quantum gravity. In 
particular we review the [P, Q] = 1 equations and combine the reasons for expecting 
the solutions to have similar non-perturbative sicknesses to the [P, Q] = 1 solutions 
of pure 2D quantum gravity. Next we derive the most general string equation 
compatible with an exact Boussinesque flow symmetry (assuming no new dimen- 
sionful parameters appear at the non-perturbative level). We study the solution for 
vanishing magnetic field in particular. The cosmological constant z — > +oo limit is 
fixed by the established genus expansion. In the z — > — oo limit we assume that all 
singularity-free solutions have an asymptotic expansion in which case there are two 
possibilities. One of these leads to the problematic [P, Q] = 1 solution, while the 
other has p — > as previously [2]. There are at most a discrete number of solutions 

3 



with the latter asymptotic, and a numerical study reveals only one. It is real and free 
of singularities. Since our Lee- Yang (5,2) solution has the same leading asymptotics 
we display it too, for comparison. 

In section 4 we construct the most general string equations compatible with 
the qth KdV flows. As mentioned above we generalise the introduction of the a 
parameter, utilise the close relationship to the generalized hierarchies, and in sect. 
5 discuss the modifications of the Dyson- Schwinger VF-algebra constraints. 



2. Review 

In this section we review the [P, Q] = Q formulation of the (2m — 1,2) models which 
was developed in [1] [2] [3] , including a complete discussion of the role played by the 
non-perturbative parameter a, the boundary cosmological constant. Central to the 
discussion is the requirement of scaling in the models and the principle that the 
KdV flows are preserved. We derive the most general string equations compatible 
with these requirements, using the Lax pair representation of the KdV flows. This 
representation readily makes contact with Douglas' differential operator formulation 
of one-matrix models, and prepares the way for the generalisations presented in 
later sections. These equations have been shown to have real, pole-free solutions 
[1][2][18]. 



2.1 Matrix Models 

The original one-hermitian matrix models [6] provided an exact solution to the 
(2m — 1,2) models via the string equations for the string susceptibility, p, together 
with the KdV flows. The solutions for p obtained from these string equations, 
although well defined in perturbation theory, produce problematic non-perturbative 
solutions for the m-even models: The only relevant real solutions to the string 
equations possess poles. The physical interpretation of these poles is unclear, and 
their presence violates the Dyson-Schwinger equations of the models. The relevant 
pole-free solutions to the m-even equations are the triply truncated solution of 
Boutroux and its generalisations [8] [14] which are complex and therefore physically 
unacceptable. 
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The problems of the definition may be traced to an instability of the one- 
hermitian matrix model at its m-even critical points, as a careful study of the 
associated scaled eigenvalue problem reveals [10]. A complementary study of the 
asymptotic behaviour of the m-even string equations reveals the presence of real 
'instanton' solutions in the single-well eigenvalue problem. The presence of the 
instantons is not by itself a signature of instability. However, the local topology 
of the eigenvalue space and the form of the effective potential for the scaled 
eigenvalues demonstrates that the definition of the m-even critical points is unstable 
to eigenvalues tunneling into a different configuration. 

Later, by studying one-complex matrix models[17] [1] an alternative exact 
solution to the (2m — 1,2) models was constructed. A different set of string 
equations which have the same perturbation theory as the previous definition was 
found. These equations possess real pole-free solutions however, thereby providing 
a more satisfactory non-perturbative definition. The KdV flows of the earlier 
definition assume a central role in these models, since it turns out that the equations 
are the most general consistent with this structure. 

The stability of these models may be traced back to the local topology of 
the scaled eigenvalues: The one-complex matrix models studied were formulated in 
terms of the combination MX f M and the new solutions had positive definite scaled 
eigenvalues. The local topology of the critical theory is thus R + in contrast to the 
1R of the original one-hermitian matrix models. The resulting 'infinite wall' at the 
boundary has the effect of removing the eigenvalue tunneling problem. The same 
effect may be obtained simply by imposing an M + topology on scaled eigenvalue 
space. 

The differential of the string equations of the original one-hermitian matrix 
model definition of the (2m — 1,2) models may be written as the canonical 
commutation relation between the operator Q representing position X s and the 
operator P representing momentum [P, Q] = 1. The operators P and Q are 
differential operators[7] in the scaled parameter z. A3 In particular, Q = d\2 + w 2 
where ui = —p and d = d/dz. 

A3 For unitary models, z is the cosmological constant and hence couples to the 
puncture operator V. 
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For a one-matrix model defined on a half-line (i.e. the scaled eigenvalue 
space has topology R+) we may write the differentiated string equations as the 
canonical commutation relation between Q and the relevant conjugate momentum, 
P, which is scale transformations A s ^- about the wall at A s = 0: [P, Q] = Q. This 
represents a sort of gauge fixed version. For full generality we should introduce 
another parameter into the theory: the scaled position, cr, of the wall. Thus the 
canonical momentum is (A s — cr)-^- = P — aP. We thus have: 

[P-aP,Q] = Q-a (2.1) 

which represents the differentiated string equations of the stable definition of the 
(2m — f , 2) models. 



2.2 The [P, Q]=Q Definition of the (2m - f , 2) Models 

In Douglas' differential operator prescription for the (2m — 1,2) series, local 
operators in the theory are constructed via fractional powers of Q, the coordinate 
operator A4, giving the infinite set Ok ~ Q\k + | , k = 1, 2, . . . and their dimen- 
sionful couplings Operator insertions are structured according to the KdV flows: 

^ = K [Q + \k + ^,Q] (2.2) 

Here to an d the cosmological constant z are seen to be related by the non-universal 
normalisation k by setting k to zero in the above: Kto = z. Thus in the unitary model 
((3, 2) = pure gravity) where z is the cosmological constant, we have QX^ + ~ V the 
puncture operator. In what follows, we shall normalise the KdV flows with k = —1. 

The string equations realising the commutation relations (2.1) may be derived 
using the principle that the KdV Hows are preserved even beyond perturbation 
theory. 

The scale transformation operator P is simply 

P=J2 Xoo(k+ht k ^- (2.3) 



A4 See section 4 for a brief review of the fractional powers of differential operators. 
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where, setting the scaling dimension of Q to 1, the dimensions of the t k may be 
derived from (2.2). Using (2.3) and (2.2), the translation operator follows as 

1 d 
P=Y / Xoc(k + -)t k - 

k=l 2 dt k~l 

Using these definitions and (2.1) the differentiated string equation for the (2m — 1, 2) 
models is: 

[]T Xoo(k + ^)t k Q+Xk + ^-a £ Xoo(k + ^)t k Q+Xk-^-^d,Q]=Q-a (2.4) 

fc=l " ~ k=l 

where Q = dX2 + u 2 . Using the fact that [Q + Xk + |, Q] = 1Z\'[— u 2 ] k +i equation 
(2.4) is a differential equation for u 2 : 

-K\"' + (u 2 - a)K\' + \u\' 2 K = (2.5) 

where 

TZ = Xoo(k + -)t k TZ k - z 
k=i 1 

The 7^fc's are the Gel'fand-Dikii differential polynomials A5 in — u 2 - In the above, we 
have used the recursion relation ViR, k +i = V 2 lZ k where V\ = d and V 2 = |<iA3 + 
u 2 d+ ^u\' 2 . The requirement for them to vanish at u 2 = fixes them uniquely up 
to the normalisation TZq, which we set to 2. When multiplied by TZ (2.5) may be 
once integrated to give: 

(u 2 - a)TZX2 + -KKX" - hllX')X2 = (2.6) 

where the matrix model tells us to fix the constant of integration to zero by the 
requirement that in the z — > +oo limit we must have the asymptotic expansion of 
7Z = coinciding with the hermitian matrix model perturbative physics. It is easily 
verified that the string equation (2.6) has the generalised Galilean transformations 

u 2 — > u 2 + e 

a — > a + e 

z^z + e^tt ( 2 ' 7 ^ 
3 

t k ^t k - e(k + -)t k+1 k>0 



A5 In the original work [21] the differential polynomials i?fe[w] are normalised such 
that 4R k = TZ k . 



as a symmetry. Using this symmetry we may perform a redefinition of the t k s 
in order to set a to zero. This corresponds to putting the potential wall in the 
matrix model at the origin of eigenvalue space. The resulting string equation was 
discussed in refs.[l] and [2] where it was argued that it possesses real, pole-free 
solutions for the mth model with the asymptotics — > zXl/m (0) in the z — > +00 
(—00) limits. A numerical solution for pure gravity (m = 2) is presented in [2]. 
Further analytical and numerical study has demonstrated the consistency of KdV 
differential flows between all the m-critical models, and the m = 1 and m = 3 
solutions were displayed [18]. The symmetry (2.7) may be rewritten as a "flow" for 
U2 under the parameter a: 

£ = i+£*»(*+!>^ = -™' < 2 - 8 > 

where we have used the KdV flows for v>2- dt k U2 = —TZX't+i in the last step. 
Equation (2.8) is consistent with the interpretation of the differentiated string 
equation (2.4) or (2.5) as a scale-invariance equation for n?,'- 

v-* * /, In du2 z 8U2 du2 „ , ^ 

gAoo( fc + -) it — + -— + a— +ua = (2.9) 

Equations (2.8) and (2.9) may be written as constraints on the partition function 
of the theory, which plays the role of a tau function, r, of the KdV hierarchyA6. 
They are then seen to be the familiar L_i and Lq Virasoro constraints, but with 
modifications in the presence of arbitrary a: 



(2.10) 



L_it = ]T Aoo (k + - )t k — + -tX2 r = — 

j^i 2 dtk-i 4 da 

> /r 1, dr 1 dr 

L°T=Z\00(k+-)t k — + -T = *- 

The rest of the Virasoro constraints are similarly modified: 

> /, In, dr 1 > 9X2t , dr 

L nr= 2^ Xoo(k + -)t k — h-^ART: w: = crAn+l— n>l 

fc =o 2 dt k+n 4 ^ dt k -!dt n - k da 

(2.11) 



A6 r and U2 are related by U2 = 2dX2 In r 
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That a produces a term of the form aXn + ld/da in the L n is consistent with 
the fact that the Virasoro constraints represent diffeomorphisms in the space of 
eigenvalues of the one-matrix model. The terms in a then arise naturally as 
boundary terms from infinitesimal variations of the position of the "wall" at A s = a. 

It was also noted in ref.[3] that a plays the role of a boundary cosmological 
constant in the theory. It takes on the mantle of the combination £ m _i/(m + 
^)t m which was identified as such for the m-critical [P, Q] = 1 theory to first order 
in tk, k < m[19]. Here the identification is not perturbative. It is also cleaner 
and more natural since a couples directly to which through (2.7) is seen to 

lead to an eX—a£ dependence for the macroscopic loop u(£) ~< e\£(u2 + d\2) >, 
identifying a as the boundary cosmological constant. Thus L_i is the conjugate 
i.e. boundary length operator — as is already obvious in the one-hermitian matrix 
model from the corresponding Ward identity [19]. We are taking consistently the 
one-hermitian matrix model definition for a macroscopic loop; these issues and 
the full loop equations are discussed in depth in ref.[4]. The above observations 
lead naturally to generalisations when we consider the general (p, q) model. These 
are discussed, together with a derivation of (2.11) and (2.10) through the Dyson- 
Schwinger equations, in section 5. 

3. The Ising model 

In preparation for the generalisation of the [P, Q] = Q formulation to all of the 
(p,q) models, we study the Ising model (4,3). We derive the most general string 
equation compatible with the Boussinesque flows, but postpone the now subtle issue 
of boundary parameters until sections 4 and 5. We begin with a brief review of the 
[P, Q] = 1 formulation of the Ising model [20]. 

3.1 [P,Q] =1 Ising Model 

Consider the following two-hermitian matrix partition function: 




(3.1) 
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where 



S(M+, M_, H, g, c) =Tr(M + A2 + M_A2-2cM + M_ - geXHM + XA- geX-HM_XA)l 

(3.2) 

This defines the Ising model on a random surface where M± represent the two Ising 
spin states at the vertices of the diagrams in the usual large- N expansion of (3.1). 
H is the magnetic field. 

Following refs.[15] and [16] the partition function may be expressed in terms of 
the norms h n of certain polynomials V n X± which are orthogonal with respect to a 
measure weighted by the potential (3.2): 

JdXdp, eX-^S(X, /j)V n X + (X)V m X - (/j) = h n 5 nrn (3.3) 

The V n X± satisfy a recursion relation: 

XV n X ± (A) = V n+1 X ± (A) + R n X ± P n _xA ± (A) + S n X ± V n - 3 X ± (A) (3.4) 

The string equations arise as the double scaling limit [20] of identities derived using 
(3.4) and (3.1). The scaling functions U2 and us (related to the string susceptibility 
and the magnetisation), together with the variables u, \i and B (the physical string 
coupling, cosmological constant and magnetic field) arise as the scaling parts of the 
quantities in the recursion relation. They are related to the free energy V as follows: 
U2 = — 3/2z/A29A2 At r and us = z/A2<9 m <9b]T '. Henceforth we shall absorb v into the 
quantities fx and B defining z = fi/v and B = B jv. 

Another approach to derive the string equations is to study the double scaled 
limits, P and Q, of operators P mn and Q mn defined by: 

Qmn= dXdiueX-^-S(X,iu) -= \i -= 

J 1 Vh m Vh n 

and 

P h\J \ N c( \ \ ]jrnX + (A) d V n X-(jl) 

Pmn= dXd/j,eX--S{X,n) -= -= 

J 7 Vh m dn ^/h n 

The string equations then arise from the requirement that P and Q, which are 
differential operators in z, satisfy [P, Q] = 1. This approach was first proposed in 
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[7] and carried out explicitly for the two-matrix model in [22] to yield the string 
equations for the Ising model. A7 

In the double scaled limit, the operator Q becomes a third order differential 
operator in z: 

Q = d\3 + ^{«2, d} + u 3 

Here, d denotes d/dz. The critical point for the Ising model is realised when P is a 
fourth order differential operator. The requirement that P and Q satisfy [P, Q] = 1 
fixes P to be Q + \4/3 where the '+' denotes the differential operator part of the 
pseudo-differentialA8 operator QA4/3. 

A simple calculation yields the following string equations: 



^(«A"" 2 + Vu 2 u\" 2 + ^ 2 A')A2 + 6u 2 X3 - 8u 3 X2) = z 

\{uX" 3 + ?>u 2 u 3 )=B 



(3.5) 



where we have integrated once with respect to z. One integration constant has been 
absorbed into z, and the other is identified with the magnetic field B [7]. (To match 
some of the conventions used in refs.[20] we must redefine z — > z/2, v — > u/y/Q and 
u 3 — > — u 3 and recall that the string susceptibility p is —u 2 -) We may study some of 
the tree level physics obtained from these equations by taking the leading behaviour 
for large cosmological constant z. This amounts to neglecting all derivatives: 

1 2 

— -uA3 2 H — -uA23 = z , 

2 3 (3.6) 

2u 2 u 3 = B 

Eliminating u 3 from the resulting equations, and adopting the above conventions, 
we have the following equation for the string susceptibility at the sphere level: 

pX3+ 3^2 =Z (3 - 7) 



A7 The authors in ref.[22] also studied potentials of higher order than quartic to 
discover that the two-matrix model can be tuned to critical points other than the 
(*, 3) models. This issue will not concern us at present. 

A8 These objects are briefly reviewed in section 4. The reader is referred to one of 
many fine works on the calculus of pseudo-differential operators such as ref. [23] [24] 
for a comprehensive treatment of the subject. 
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From this we may derive expressions for tree level correlators of the puncture 
operator V and the spin operator S. 

The Ising model (3.1) with H ^ is not symmetric under the interchange of 
M + and M_ . Such an interchange translates into d — > — d in the differential operator 
formalism[20] and this is the origin of the Z 2 -odd transformation properties of B 
and us. Studying the Zi symmetric sector by setting B and 113 to zero, we have the 
string equation for the string susceptibility: 

jL p X"" - i( P A')A2 - p P \" + pX?> = z (3.8) 

Large z expansion supplies the genus perturbation theory which is uniquely deter- 
mined from the the spherical contribution p(z) = zXl/3. 

Equation (3.8) is very similar to the string equation for the Lee- Yang singular- 
ity (the (5, 2) model [25]) where instead the coefficient of p\"" is 1/10. A numerical 
solution to that equation with the asymptotics p — > ±zXl/3 for z — > ±00 was 
presented in [26]. The family of (2m — 1, 2) string equations was studied in ref.[14] 
in order to construct the solutions analogous to the 'triply truncated' solutions 
found by Boutroux for Painleve 7 . In that work, the authors demonstrated that the 
generalised truncated solutions were real for the m-odd models and complex for the 
m-even models. The Lee- Yang string equation falls into the former category and 
the Painleve 1 equation into the latter. 

It is a simple matter to repeat the analysis for the Ising equation (3.8). The 
reality of the truncated solutions is determined by the nature of the solutions to the 
associated polynomialA9 sX2 — As + 3A = 0, where 1/A = 2/27 for the Ising model 
and 1/10 for Lee- Yang. If s is complex then the truncated solutions will be real, 
and for real s they are complex. We see that s is complex for A > 1/12 and so the 
Lee- Yang solution is real (s = 5 ± iy/E) and the Ising model solution (s = 27, 9/2) 
is complex. 

The alert reader will note that the same polynomial determines the reality 
of the 'instantons' in the asymptotic expansion of (3.8) when addressing the 
question of Borel resummability [12] [13]. These 'instanton' solutions are merely the 
leading exponential corrections to the perturbation expansion as a representation 

A9 We refer the reader to ref.[14] for details. 
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of the full non-perturbative solution. The real instantons and corresponding non- 
resummability for the Ising model is purely a consequence of the unitarity of the 
theory (whereby all terms in the genus expansion contribute with positive sign) and 
the typical (2n)\ growth of the perturbative series. 

Later we will find that the [P, Q] = Q formulation supplies the same genus 
expansion as for (3.8) and hence the same resummation properties for positive 
cosmological constant. Nevertheless, we shall explicitly have a real non-perturbative 
solution. This situation was already discussed for the (2m — 1, 2) models in [2]. 

3.2 [P, Q]=Q Ising Model. 

We must first recall the underlying structure which exists perturbatively for the 

[P, Q] = 1 definition of the (*, 3) models, the Boussinesque hierarchy, which defines 

the flows of Q = dX3 + (3/4){u 2 , d} + u 3 : 

BO 3/r -+- / 

< =iAQ\-^i^ ,Q] 1 = 1,2; fc = 0,l...oo (3.9) 



dt^ k '"' 3 + 

The ti^k are an infinite set of parameters which parametrize the hamiltonian flows 
of Q. k is a non-universal normalisation parameter. We may write the equation 
(3.9) as a pair of equations for u 2 and u 3 : 

a (t) = KVXijxKXji^+i = KVXij 2 'JlXji t k i, 3 = 2, 3 (3.10) 

The lZXji ; k are differential polynomials in u 2 and u 3 , and ct( 2 ) = 3/2, = 1. The 
subscript bracketed (i) indicates "no sum on i" . The TlXji^ are the generalisation of 
the Gel'fand-Dikii differential polynomials for the KdV hierarchy. The equivalence 
of the two hamiltonian structures defined in (3.10) implies a recurrence relation 
between them. Requiring them to vanish at u 2 = u 3 = fixes them completely up 
to the normalisations IZXii^. The first few are: 

ftA2i i0 = 3; ftA3i i0 = 0; 

TZX2 2)0 = 0; ftA3 2i0 = 3; 

3 

ftA2 M = 2w 3 ; ftA3i,i = -u 2 ; 

TZX2 2)1 = -^(uX" 2 + 3w 2 A2); ^A3 2 ,i = 2u 3 ; 

TlX2 h2 = -^uX"" 2 - ^u 2 uX" 2 - ^(u 2 X')X2 - ^uX3 2 + ^uX2 3 ; TZX3 h2 = ^(uX" 3 + 3u 3 u 2 ) 
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The objects VXij 1 and V\ij 2 define the first and second hamiltonian structures 
of the Boussinesque hierarchy. They are a shorthand notation for the fundamental 
structures defining the Poisson bracket for functionals of u 2 and u 3 . The explicit 
expressions for them are: 

2 1 
£>A22 2 = -d\3 + -u 2 \ + u 2 d 

2 

£>A23 2 = u 3 d + -u 3 X' 

VX32 2 = u 3 d+-u 3 X' 
3 

X?A33 2 = ~^d\5 - —u 2 dX3 - §w 2 A'dA2 + (-Jw 2 A2 - |w 2 A")d + (--u 2 u 2 X' - ^u 2 X'") 
lo 12 o 2 o 2 12 

and £>A22i = £>A33i = 0; £>A23i = £>A32i = d. From the equation 

dxi ' 
crti,o 

we make the identification z = nt\fi. In what follows we shall set k = —1. The 
scaling dimension of Q supplies a natural length scale in the theory. Fixing its 
dimension fixes the scaling of the t\^. If we assign the scaling dimension 2 to u 2 , 
we obtain [Q] = 3, [t^k] = — (3k + /), and [u 3 ] = 3. 

The central assumption which leads uniquely to the string equation is that 
the Boussinesque Hows (3.9) hold at the non-perturbative level. We construct P, 
the generator of scale transformations in the theory, out of the parameters in the 
Ising model with magnetic field: ti )2 (which defines the (4,3) model), £1,0 = —z, and 
£2,0 oc B. 

^ d d „ d 

Using [d tl 2 , Q] = — [Q+A7/3, Q], we have for our differentiated string equation: 

-3Q = 



-7 tl , 2 Q + \7/3-z^-2BjL,Q 



which is a pair of scaling equations for u 2 and u 3 : 

(dui dui ^dui \ 
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Identifying B = — 2t 2 ,o? these equations may be succinctly written as: 



V\ij 2 Tl\j = 



(3.11) 



where 1ZX2 = TZX2i j2 — z and 1ZX3 = 1ZX3i^ — B and we have set ti^ = 1/7. It 
should be noted here that we may write the derivative of the [P, Q] = 1 string 
equations (3.5) as 



Finally, to obtain the string equation, we multiply equation (3.11) on the left by 
TZXi giving us a total differential and integrate once with respect to z, to give: 



(For convenience of notation we have exchanged the superscripts on the lZX€s for 
subscripts.) We have set the constant of integration to zero by requiring that our 
perturbative physics obtained in the z — > +oo limit is the same as that obtained 
from the matrix model via equations (3.12). Indeed, with the constant in place, the 
tree level string equation is: 



from which we may obtain the equations (3.6) by setting each bracket and the 
constant to zero. We can then obtain the same tree level physics as the [P, Q] = 1 
definition. With the constant set to zero we must always follow this procedure at 
any level of perturbation theory to match the physics of the [P, Q] = 1 equations. 
This is a direct consequence of the fact that the structure of the string equation 
(3.13) always admits a 7Z2 = TZs = solution. 

3.3 The Solution for the String Susceptibility. 

We now study the physics of equation (3.13) in the absence of the Z2 breaking 
quantities 113 and B. Setting them to zero and adopting the conventions of (3.8) we 
have the following string equation for the string susceptibility, p: 



VXij^RXj = 



(3.12) 




( — U2X3 H — uX2 3 — z)X2 — 112(2^2^3 — B) = constant 




(3.14) 
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where 



Uj = ^pX"" - ^( P X')X2 - ppX" + pX?,-z 



We now have an equation for p which is structurally identical to that obtained for 
the string susceptibility of the Lee- Yang model (5, 2): The coefficient of the first 
term in (3.14) would instead be 1, and the expression for IZly would have 1/10 as 
the coefficient of the first term. The two equations are similarly analysed. Using 
dimensional arguments the asymptotic expansion for z — > +oo may be shown to be 
of the form 



where the A^s are dimensionless constants. By substitution these can be seen to 
be determined uniquely once A is known. Since the same is true of the [P, Q] = 1 
string equation (3.8) it follows that the resulting perturbative expansion is identical 
in both cases once we have set A = l. As in ref.[2] we assume that all real solutions 
without an asymptotic expansion in the z — > — oo limit have poles. Requiring an 
asymptotic expansion in this limit we find that the sphere solution is either p = 
— \z\Xl/3 or p = 0. It is clear from the above discussion that the first choice leads 
to the same perturbative physics as the [P, Q] = 1 definition and probably the same 
problematic non-perturbative solution. (Certainly there are at most a discrete 
number of solutions in this case with the latter being one of them A 10.) This 
expectation is reinforced by the matrix model understanding of the (2m — 1,2) 
models; in particular the Lee- Yang model's [P, Q] — Q definition has the spherical 
physics fixed in the one-matrix model [1]: The end-point of the eigenvalue density 
of the model corresponds to the string susceptibility p in the spherical limit. In the 
z — > +oo limit the endpoint of the density pulls away from the wall where we recover 
locally the 1R topology and hence the physics is identical to that of the [P, Q] = 1 
models; po = zXl/3. In the z — > — oo limit, the endpoint pushes up against the wall, 
and po vanishes. 

The appropriateness of the choice p = in the z — > — oo limit for the string 
equation (3.14) is made manifest by a detailed analytic study of the family of 

A10 A numerical study of the string equation, using the techniques described later, 
failed to find a pole-free solution with these asymptotics. 



p = zX- ^XooA l 
6 i=0 




Xi 
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equations of this form and their solutions [18]. The KdV and Boussinesque flow 
structure of the (2m — 1,2) and (*, 3) models respectively, may be shown to preserve 
the monodromy data of the linear problem associated to the string equations if and 
only if the po — asymptotic is preserved. This asymptote is fixed in pure gravity [2]. 
Note that (2,3) = (3,2), a fact which is trivially verified using the expression for 



Using the same procedure as for the z — > +oo limit, in the z — > — oo limit with 
Po = we have the asymptotic series 



The -B/s are again fixed uniquely by the initial Bo = —4/9. 

We continue by studying linear perturbations e about the leading behaviour 
of the string susceptibility in the large z regions. In the z — > +oo limit we have 
po = zXl/3 and so we try to find a solution for e(z) by substituting p = zXl/3 + e(z) 
into equation (3.14). Following the WKB prescription for large z we expect the 
exponential behaviour e~eA— f(z) with |/| >> |/A'| >> |/A"| ... so that eA(n)/e^ 
(-fX')Xn. Using this and keeping only leading order we find that 



and lZX(n)j(po + e) rs (-fX')XnlZi(po + e). (Derivatives of p do not survive in this 
limit.) The string equation (3.14) becomes 



From this we find e(z) = AeX-6/7azX7 /6 with aA2 = 9/2, 9 or 9/4. As three of these 
solutions are exponentially growing perturbations we must set their coefficients to 
zero. 

In the z — > — oo limit we study perturbations around the leading non-vanishing 
behaviour for the string susceptibility, the torus term p = —A/9zX2. This time, as 
po and it derivatives are subleading, we have 1Zi(po + e) ~ e2/27(/A / )A4 — z. The 
string equation then becomes, to leading order 9/16zA2 + z/27(fX')X6 = which 
gives (/A')A6 = — 243/16|z|. This yields the following solution for the exponential 



TZX2 2A . 





(A( /A ')A4 - (fX')X2zX± + 3zX^j X2 (^zX± - (fX')X2^j = 
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corrections in the z -> -oo limit: e(z) = BeX-6/7(3zX7 /6 with (3X6 = -243/16. 
Again, three of these solutions have positive real part and their coefficients must 
be set to zero to match the chosen asymptotes. Six integration constants have now 
been determined locally in a sixth order differential equation and so we expect at 
most a discrete number of solutions with the above asymptotics. 

Further progress was made with numerical techniques identical to those used in 
[2] to find the solution for the [P, Q] — Q pure gravity model (m = 2) and in [18] for 
the m = 1 and m = 3 models. The solution of the differential equation was treated 
as a two-point boundary value problem. We employed a NAG FORTRAN library 
routine (D02RAF) to solve the problem by relaxation. The program uses Newton 
iteration with deferred correction, and allows user specification of the initial mesh 
and approximate solution. The absolute error tolerance was set at ~ 10A— 5. 

We chose to solve the string equation differentiated once, removing a factor 
of TZj. This allows for a more numerically well-behaved highest derivative, since 
otherwise the expression for the highest derivative obtained from the string equation 
contains factors of 1/lZi. It was ensured that the correct solution was found by 
including more terms in the asymptotic series to calculate the boundary conditions. 
The stability of the solution was tested by performing the integration with a number 
of different values for the boundary, e.g. z = ±100 and z = ±10. 

The non-perturbative solution for the string susceptibility of the Ising model 
is displayed in figure 1(a), where the integration was carried out on a mesh of 
1600 points in the range ±200. For contrast, in figure 1(b), we display the non- 
perturbative solution for the Lee- Yang model. 



4. The (p, q) String Equations and the qth KdV Hierarchy. 

We present here the generalisation of the [P, Q] = Q formulation to the {p, q) 
minimal models, deriving the unique string equation consistent with the requirement 
of preservation of the qth. KdV flows. A parameter analogous to a is included in 
the discussion, which in the Ising model is seen to be the boundary magnetic field. 
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4.1 The qth KdV Hierarchy. 

We begin by reviewing the basic tools of the formalism, the KP hierarchy and its 
(/-reductions, referred to as the gth KdV hierarchies. The KP hierarchy may be 
formulated in terms of the pseudo differential operator L = d + J2i=i Xoofi(t)dX—i, 
where dX—1 is defined by dX—lf = J2j=o Xoo(—l)XjfX(j)dX—j — 1. Integer powers 
of L generate a basis for the complete set of objects which commute with it. Taking 
the differential operator part of these (denoted by a '+' subscript) generates a set 
of evolution equations for L: 

^-=k[L + Xt,L] (4.1) 

The t r parametrise the infinite set of flows thus defined. Equation (4.1) defines the 
KP hierarchy. 

The qth reduction may be constructed in terms of the object Q = LXq and then 
requiring that Q- = 0: 

wr 4Q + x ? Q] <42) 

When r = mod q the flows of equation (4.2) are trivial, so we modify our notation 
explicitly to highlight the values of r which are mutually prime with q: 

()Q = k[Q+X^±^, Q] I = 1, 2, . . . q - 1; k = 0, 1, . . . oo (4.3) 



The indices / and k now span the set of non-trivial flows r = qk + 1. Equation (4.3) 
defines the qth KdV hierarchy. The differential operator Q may be written in the 
form: 

Q = dXq + Xqcti{ui, dXq — i} (4-4) 

i=2 

and (4.2) defines a pair of Hamiltonian equations for the {ui}: 

dui 



= K{Hi ! k+i,u i }i = K{Hi }k ,Ui} 2 (4.5) 



dt ljk 

The hamiltonians are constructed from fractional powers of Q in the following way: 



^ = kqTl J^Q X — d z = kqTl^- 



19 



where the residue of a pseudo-differential operator is simply the coefficient of the 
dX—1 term. We note here that the Hi,k are not well defined for the relevant solutions 
Ui, which grow as a power of z as z — > +00. Their explicit appearance here uncovers 
the structure of the q-KdV system. They themselves are never used with these 
solutions, both here and later. 

These q-KdV systems are 'bi-hamiltonian': They possess two Poisson brackets 
between functionals W, V of the {ii;}: 

r SW 5V 

{W[UI V[u]} h2 = Jtedy — iuiW^y)}^— (4.6) 

The fundamental Poisson brackets {iii(x),Uj(y)}i j2 may be written 

{iti(x), Uj(y)} li2 = V li2 \ij(x)5(x - y) (4.7) 

The objects T>Xij 1 and T>Xij 2 are a set of differential operators. Using them we 
may develop (4.5) further 

= K{Hl,k+l,U i } 1 = K{Hl,k,U i }2 



dt 



Lk 



..8Ht,k+i ^ x . .6Hi,k (4.8) 

OUj OUj 

> V 1 Xij1ZXji, k+1 = V 2 XijTlXji jk i, j = 2, 3, . . . , q 



where the IZXji^ are differential polynomials in the {ui}. They are the generali- 
sation of the Gel'fand-Dikii differential polynomials encountered in the q = 2 case. 
The last line in (4.8) is a recursion relation among them. Requiring them to vanish 
at {u^ = fixes them uniquely, up to the normalisations IZXi^o oc qSXi — 1 /. In what 
follows we set these at TZXi^Q = q5Xi — 1/, and we set the overall normalisation k 
to -I. The second Poisson bracket in (4.7) is in fact the WX(q) = WA q -i -algebra 
(at a particular value of the central charge) where the {ui} are the q — 2 currents. 
In particular, u 2 corresponds to the energy-momentum tensor. For example in the 
q = 2 case we have: 

{u 2 {z), u 2 (y)} 2 = (^dXS + u 2 d + ^u 2 X'^j S(z-y) 
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which is the WX(2) or Virasoro algebra, and for the case q = 3 we have a WX(S)- 
algebra: 

/2 1 \ 

{u 2 (z), u 2 (y)}2 = ( g^A3 + 2 U 2A ; + it 2 d) <5(z - 2/) 

{1*2(3), u 3 (y)} 2 = [u 3 d + g^A'^ <J(z - y) 

{1x3(2;), «2(2/)}2 = (u 3 d + 7^u 3 \'^ S(z - y) 
{u 3 (z),u 3 (y)} 2 = 

= ("^A5 - ^u 2 dX3 - ^u 2 X'dX2 + (-^u 2 X2 - ^u 2 X")d+ (-^u 2 u 2 X' - ^A'")) S(z - y) 
which forms the second hamiltonian structure for the Boussinesque hierarchy [27]. 

4.2 The {p, q) String Equations. 

As discussed before, the matrix model formalism motivates us to work with the 
operator Q which plays the roole of the continuum limit of a position operator 
in the orthogonal polynomial basis. We construct a realisation of the equation 
[P — crP, Q] = Q — a in terms of the parameters t^k- Heat kernels of Q generate 
the correlators of the observables in the theory, which are the macroscopic loops. Q 
therefore supplies a natural length scale in the theory. From the definition of Q in 
(4.4) and the q-KdV flows in (4.2) the scaling of the {ui} and {ti^} may be deduced. 
If we assign the scaling dimension 2 to u 2 we obtain the following dimensions: 
[ui] = i; [t^k] = —(qk + /); [Q] = q. Using these, we construct the generator of scale 
transformations, P: 



P = E a? - 1 E Xo °(i k + l )hk 



2=1 k=0 



dti jk 



and the generator of translations, P: 



P = E A? - 1 E Aoo(q/c + l)t hk 



d 



1=1 k=i 



dti^-i 



By then adopting the principle that the qth KdV hierarchy holds we have 



[E A? - 1 E Aook(q/c + l)ti, k QXk + - - aQXk H 1 




,Q]=q(Q-<j) 



+ 



(4.9) 
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We may rewrite this as a set of scaling equations for the {ui}: 

■ du ■ du ' \ 

(4.10) 

The q — 1 objects t m o are proportional to the parameters coupling to the relevant 
operators in the theory, the O m $: In particular we have = —z, 2t 2 ,o = — etc 
(after setting k = — 1), and so in rewriting (4.9) we must use the identity 

and in order to interpret it as a scaling equation we have made the following 
identification: 

a {l) -^ = -qVXij i n\j (4.11) 
The scaling equations (4.10) may be written succinctly as 

(VXij 2 - aVXij^TZXj = (4.12) 
where the objects TZXi in the above equation are 

IZXi = Y M - 1 Xoo(k + -)ti,kR<\ii,k 

1=1 fc=0 1 

= J2 X( i~ 1 J2 Xo °( k + -)ti,k^Xii^ + (i - l)*»-i,o 
z=i fc=i ^ 

In the above, we have used that IZi^Xj = q5Xj — l L . Equation (4.12) is the 
differentiated string equation. Its structure is an explicit realisation of the WX(q)- 
algebra structure inherent in the second Hamiltonian structure of the q-KdV 
hierarchy. In analogy with the case explicitly worked out from the matrix model, 
the string equation is obtained by multiplying on the left by TZXi giving a total 
derivative, and integrating once with respect to z. (The integration constant is 
then set to zero using perturbation theory. See section 3.3) 

It should be noted here that the a-deformed differentiated string equation for 
the (2m — 1, 2) models may also be written in the form of (4.12): (D 2 — aVi)TZ = 
It is now apparent that the process of introducing a into the formalism may be 
regarded as forming a linear combination of the second hamiltonian structure T>2 
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and the first, V\. This is always possible as the two structures are 'coordinated', 
in the sense of ref.[24]. In this picture (4.12) is indeed the natural generalisation of 
the equations first found for [P, Q] = Q definition of the (2m — 1,2) models. 

That 1ZXi(VXij 2 — crVXij^TZXj is always a total derivative must be demon- 
strated. The integral JdzTZXiVXij 12 TZXj is 

J2 Xq-l Xoo(k+ -)(n+—)ti tk t m ,n [fcXii tk T>Xij li2 'R,Xj rnin dz 

l,m=l k,n=0 Q Q J 

We then have using equations (4.6) and (4.7): 

JllXii,k'DXij 12 KXj rn , n dz = Jdz dyKXii, k (z)VXij 12 {z)8{z - y)H\j m ,n{y) 

= J dzdy 5^ {u ^ U ^ 2 ^) 

= {7~Ll,k, 7~Lm,n} = 

where the last line is simply the statement that all the Hamiltonians in the q-KdV 
hierarchy are in involution. However, the above equation holds for all functions Ui 
satisfying the boundary conditions (which are those neccessary for the Hamiltonians 
H.i,k to be well-defined). This can only be true if the integrand on the left hand 
side of the equation is a total ^-derivative. Thus we conclude that we may always 
integrate (4.12) to obtain the string equation. 

As an example, we have the string equation for the (*, 3) models: 

lu 2 n 2 X2 + ln 2 nx" 2 - \(nx' 2 )X2 + ( U3 - a)n 2 n 3 - 4 (n 3 nx(4) 3 - nx' 3 nx'" 3 

Zoo lo V 

- ^ (u 2 n 3 nx" 3 - ^u 2 (nx' 3 )X2 + ± U2 x'n 3 nx' 3 ^ - ^ (3« 2 A2 + u 2 x") nx2 3 = o 

(4.13) 

where 

Z=l,2fc=l 6 

and 

n 3 = S Xoo ( k + 4)*i,fc^ 3 i,fc - B 

1=1,2 fc=l 6 

(For convenience of notation we have exchanged the superscripts on the IZXVs for 
subscripts.) From this, the Ising model is obtained by setting t^k = 3/78Xli5Xk 2 . 
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5. The W— algebra Constraints. 

In this section the constraints for the [P, Q] — Q formulation of the (p, q) models are 
derived. We complete the discussion of the parameter a in the (2m — 1,2) models 
in terms of the algebra of constraints, which provides the appropriate method for 
generalisation to the (p, q) models. We discuss the possible significance of the 
analogous q — 2 extra parameters. 

5.1 The (p, q)-Model W-constraints 

In the [P, Q] = 1 formalism, the string equation leads to an infinite set of constraints 
on the partition function of the theory, which plays the role of a r-function of the 
associated gth reduction of the KP hierarchy. The [P, Q] = 1 string equation (3.12)is 
equivalent to the L_i constraint on the r-function: 

VXij-JIXj = ddXi (^j?— ) = 

(For what follows we shall work with all the £ r 's which parameterise all of the flows, 
including the trivial ones [28]. The IZXi's from the previous section are now defined 
as J2 r =i Xoo^t r 1ZXi r , and in equation (5.1) we have dXi = d/dU-i. We also have 
Ui = (i / 2) ddXi (In r).) 

In ref.[28] it was shown that the constraints derived from the string equation 
may be written as: 

WX(k) n .r = 0, n>-k + l, k = 2, 3, • • • , q (5.2) 

where the WX(k) n are the nth Fourier modes of the W X (q) -algebra generator with 
spin l k\ For example WX(2) is the stress tensor and its Fourier modes (usually 
denoted L n ) satisfy the Virasoro algebra. The constraints L n .r = n > — 1 then 
form a consistent set in the sense that no further constraints upon r are generated 
using the commutation relations of the modes. 

In this section we will show that equations (3.11), obtained from the [P, Q] = Q 
string equation by differentiation, imply the following constraints: 



WX(k) n .r = 0, n>0, k = 2,3,---,q 
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(5.3) 



which also form a consistent set. These constraints follow from the differentiated 
string equation, which is equivalent to the L$ constraint: 



V\ij 2 n\j = ddXi = 



where 

x d 

L = Xooiti- h const. 

i=i vti 

Using the techniques and notationAll of ref.[28], the constraint Lq.t = may be 
written: 

res x (Xd fI X2X(t, A, fj)^) .r = (5.4) 

From equation (2.8) of ref.[28], 

res\(XuX*(tX' , X)d\ui(t, A)) 

'X(t, u)res x (Xd^X2X(t, A, ^)| M=A ) r(t) XX * (tX', u)r(tX') 



res,. 



r(t) r(tA') 

(5.5) 

and thus equation (5.4) implies that the left hand side of (5.5) vanishes. However, 
following ref.[28] one has 

res x (XajX*(tX', X)d x u(t, A)) = (MQXl/q(t))_ 5(x - xX') 

We have shown that the [P, Q] = Q string equation implies the condition 

(MQAl/g)_ = 0. (5.6) 

Taking powers of T = MQXl/q, it is then straightforward to show that equation 
(5.6) implies the further conditions 

(MXnQXn/q)_ = 

for n = 1, 2, . . . , q — 1. Following the discussion of ref.[28], these lead to the 
constraints (5.3). 

The above equations correspond to the case where the scaled eigenvalue space 
is taken to be R + with the boundary, or 'wall', being at A s = 0. In fact there 



All Beware of the interchange of the names of the operators Q and L in ref.[28] 
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is no reason for this restriction on the coordinates and in general one may take 
the boundary position to be A s = a. The effect of this on the string equation 
and the Virasoro constraints of the one-matrix model has already been described 
in refs.[3][ll][4]. approach most suited to There it is shown that the eigenvalue 
boundary plays the role of the boundary cosmo logical constant on the worldsheet. 
In sections 4.2 and 3.2, by identifying the eigenvalue space with the position operator 
Q in Douglas' PQ formalism, we derived the effect of a on the string equation in 
the general (p, q) case, and on the Ising model in particular. In the latter case the a 
parameter will be seen, by precisely parallel arguments to refs.[19] [4], to be the (Z 2 
odd) boundary magnetic field. This method however takes into account only one 
of the flavours of eigenvalue in a multi-matrix modelA12, each of which may have 
their own boundary. 

Rather than working directly in the continuum limit the clearest picture might 
be expected to emerge from working explicitly with a q — 1-matrix model. (For the 
one-matrix model this was done in ref. [4] where it was used primarily to derive the 
effect on macroscopic loops.) Integrating over the angular modes leaves an integral 
over q — 1 coupled eigenvalues, and the obvious generalisation is to introduce q — 1 
boundary parameters giving the position of the q — 1 'walls' in scaled eigenvalue 
space, eigenvalues restricted to being greater than their In the Ising model, a critical 
point in the two-matrix model, linear combinations of the two flavours of loop give 
the Z 2 even boundary length and Z 2 odd boundary magnetization[19]. Given the 
results of the one-matrix model, it is then natural to conjecture that here the wall 
parameters will in a similar way provide the conjugate parameters: the boundary 
magnetic field and the boundary cosmo logical constant. Note that the latter is not 
apparent in the [P, Q] = 1 KP description[19]. It is therefore an important question 
to determine whether such a parameter exists in our formulation. We will below 
identify a Z 2 even parameter, conjugate to a redundant Z 2 even operator (namely 
W- 2 ) which we therefore suggest is (perhaps non-linearly) related to the boundary 
cosmo logical constant. Curiously, we will also uncover a further redundant Z 2 odd 
parameter (conjugate to W-±) for which we have as yet no physical interpretation. 

The reason we cannot be more definite in our identifications is because we 

A12 Strictly speaking it is a linear combination[19]. 

26 



were unable to carry through a direct analysis of the two-matrix model, due to 
certain technical difficulties in the orthogonal polynomial approach (outlined below). 
Other technical difficulties nullify the standard steepest descents approach as an 
alternative method for any multi-matrix model. 

Let us now turn briefly to the two-matrix model inserting eigenvalue boundaries 
at o"i = o"2 = 0. The orthogonal polynomials are now normalised as 

/ XodXdfx eX-^-S(X,fi)V n X + (X)VX- rn (ii) = h n 5 nrn 
Jo 7 

The problem with the method arises when we consider the generalisation of 
eqn.(3.4). The reason the recursion relation (3.4) finishes at V n -3 is because the 
potential S is of order 4, as can be seen by noting that XVX + n (X)VX — n -k (p)/hn-k 
vanishes for any k > 3 on integration with the measure implied in eqn.(3.3). The 
proof follows by converting A into jnZ^I ex P(~ ^~&) + corrections, integrating by 
parts and then using (3.3) together with the definition VX — m (fi) = /iXm+ lower 
powers. Unfortunately in our case we pick up a boundary term 

^-VX - n . k (0) J q XodX eA-^S(A, 0) V n X + (A) 

for any odd k > which is therefore equal to the coefficient VX+ n -k in an infinite 
orderA13 generalisation of (3.4). The coefficients (in particular the h n ) can now be 
determined, in principle, from an infinite set of simultaneous recurrence relations 
generalising the usual construction. We conclude that the orthogonal polynomial 
approach is at best inappropriate for analysing the continuum limit. Of course these 
arguments are unaffected by choosing general positions o\, 02 for the walls or using 
a non-even potential. 

We now turn to an approach based on the continuum Dyson-Schwinger equa- 
tions. We first review the results of the one-matrix model using this approach and 
then generalise to general (p, q) . Our starting point in the one-matrix model is the 
algebra of constraints [2]: 

L n r = n>0 (5.7) 

This differs from the hermitian matrix model in that the L_i constraint is missing. 
The KdV flows on the other hand, upon which our formulation is based, are invariant 

A13 Recall that we are interested in n ~ iV — > oo for the continuum limit. 
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under transformations generated by the full set L n : n > — 1. Thus we have a 
situation reminiscent of spontaneous symmetry breaking: the 'dynamics' i.e. the 
KdV flows are invariant under the full group (generated by the L n : n > — 1) whereas 
the 'vacuum' r is invariant only under the little group generated by L n : n > 0. The 
'Goldstone boson' is a which in (5.7) has been gauge fixed to zero. Indeed there 
are now an infinite number of vacua r(a) connected by the broken generator: 

T (a)=e\aL- 1 r (5.8) 

where we identify r(0) with the r function in (5.7). The constraints on r(a) are an 
inner automorphism of those of (5.7) plus a constraint arising from a independence 
of r: 

where L' n 
and d' a 

Taking linear combinations of these constraints: 

LXa-i = -d' a 
L\o~o = L' + aLXa-i 

(5.10) 

LXai = L[ + 2aL' + a\2L\a- 1 



n>0 


(5.9) 



= eAo"L_i L n eX—aL_i 

- eXaL-i eX—aL-i . 
da 



we get LXa n = L n — aXn + 1-^ for n > —1. These corrections can be computed 
most simply as follows: 

j^t(*) = L- lT (a) (5.11) 

L n r(a) = eAaL_x (— aadL_i)An + lL n r(0) n>0 

= aXn + lL_!eAaL_! r(0) 
d 

= aXn+l—r(a) (5.12) 

OCT 

where use is made of (5.8), the notation (adX)Y — [X, Y] (for two operators), and 
the identity eX-X YeXX = eA-adX Y. 
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Our string equation (2.6) follows from integrating the modified L constraint. 
It is clear that the system is now invariant under translations generated by L_i. 
Indeed, using (5.11), the explicit formula L_i = J2k=i Xoo(k + §)*fc^fe— 1 + z\2/4, 
and differentiating, one obtains 

- }^ Xoo(k + ±)t k - 



da ^ 2 dt k -i 

which, on multiplying by an infinitesimal e, may be interpreted as the invariance 
(2.7). Clearly under a finite translation by —a we return to our original equations. 
Thus the boundary length L_i and boundary cosmological constant a are redundant 
just as was the case in ref.[19], however the symmetries and physical significance of 
our formulation are far more transparent with 

Now consider 2D gravity coupled to a general (p, q) minimal model. As shown 
at the beginning of this section, the WX(k) a constraints, where a is restricted 
to —k + l<a<0 (and 2 < k < q), are missing as compared to the hermitian 
matrix model formulation. These, together with our constraints (5.2), generate 
the full set of symmetries of the generalised KdV hierarchy. Thus we now have 
spontaneous breaking of Y,k=2 Xq(k — 1) = \q(q — 1) symmetries to which we may 
associate \q(q — 1) new parameters: aX(k) a . These parameters (and associated 
operators) will be redundant since there will be analogous symmetries to (2.7) 
which will 'gauge' them away, however we again expect that they will have physical 
significance. Generalising (5.8) we have 

r(aX(k) a ) = S(aX(k) a ,WX(k) a )r, with 5(0; WX(k) a ) = 1 . (5.13) 

Since the VF-algebra of constraints is no longer a Lie algebra something more general 
for S than exp{][] fc a aX(k) a WX(k) a } may be more appropriate. The transforma- 
tion (5.13) induces a similarity transformation on the constraints C = SCSX—1 
generalising (5.9), where C runs over the WX(k) n with n > and the new constraints 
d/daX(k) a (which are trivially zero on r). The latter constraints give rise to the 
generalised symmetries by similar derivations to (2.7); Except for L_! these are 
non-local. Since the similarity transformation preserves the VF-algebra we see that 
the W'X{k) n still form a VF-algebra while each d' aX ^ commutes with all the other 
constraints. It can be shown that, by taking linear combinations of these constraints 
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(generalising (5.10)), one can form a new simpler set of constraints of the form 
(WXa)X(k) n = WX(k) n + corrections where now n > 1 — k and the corrections do 
not involve any WX(k) m with m > 0. It may also be shown that these corrections 
can be written as a finite sum of the form J2k,a ^Wa^M^Oa where each c;X(k) a 
is a power series expansion in all the aX(k) oi with operator valued coefficients. This 
follows because the original TU-algebra with n > 1 — k forms a closed algebra of 
constraints under commutation. By construction the commutator of two of these 
WXa constraints also closes, however, as we will show by example below, these 
constraints do not form a W algebra. 

In particular, consider the Ising model (or more generally coprime p > q = 3). 
In this case we have broken 'symmetries' W-± and W-2- The most general 

case generalising (5.8) will involve some combination of these three generators 
together with three parameters. The corresponding formulas to (5.11) will be more 
complicated involving mixed combinations of parameters and operators, because 
these operators do not commute. However, we can classify the operators under the 
Z 2 symmetry which flips spin (exchanges M + with M_). Since this corresponds 
to d — > — d (cf. comments below equation (3.5)), we have QA1/3 — > — QX1/3 and 
using (4.1), it follows that t r is Z 2 -odd(even) if r is odd(even). Using the explicit 
formulas for W n and L n (see e.g. ref.[28]) we find that these are Z 2 -odd(even) 
if n is odd(even). Thus we propose that W-2 is (perhaps non-linearly) related to 
the boundary length. We do not know what role is played by the Z 2 odd W-\. As 
mentioned above, L_i is the boundary magnetization. We will now confirm this. 

Avoiding the complications of mixing let us first consider L_i in isolation. In 
this case eqn.(5.8) again applies. A direct calculation of the modified constraintsA14, 
via the method of (5.12), gives: 

WXa n = W n -(n + 2)aXn + lW- 1 + (n + l)aXn + 2W- 2 n>0 (5.14) 

with the LXa n , n > —1, as before. It is clear (by counting) from the discussion below 
(5.13) that there are no WXa-i or WXa -2 constraints; Indeed it is amusing to note 
that this is already incorporated in (5.14) if we take n > —2. This is preserved by 
the modified algebra of constraints which by explicit computation is found to be: 

A14 We use the conventions of ref.[28]. 
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[LXa n , LXa m ] = (n- m)LXa n+m 

[LXa n , WXa m ] = (2n — m)WXa n+m — (2n + l)(m + 2)aXm + lWXa n -i 

+ 2(n + l)(m+ l)aXm + 2WXa n - 2 
[WXa n , WXa m ] = --(n - m){(nX2 + Anm + mX2) + 9(n + m) + 14}LXa n+m 
+ 2(n — m) (U Xa n+m + aXn + m + 3U Xas) 
+ + l)(m + 2)(m + l)maXn + 2LAa m _ 2 

- -(n + 2)(m + 3)(m + 2)(m + l)aAn + lLAa m _i 
+ 2(n + 2)(m + l)crAn + lt/"Acr m -i - 2(m + 2)(n + l)crAn + 2UXa. 
— (n^m)} 



where UXa n may be taken to be J2k<-2 LX<JkLXa n -k + J2k>-i LXa n -kLXak- The 
corresponding (galilean) symmetry follows, cf. (2.7), from equation (4.11): 



u 3 -> u 3 + e 
-U2 — >• «2 
a — > cr + e 



2 — > z + e—ti i 
3 ' 

B^B + e^t 2 ,i 
U,k -»■ - e(fc + 1 + ^)*i,fc+i 



fc>0;/ = l,2 



and the string equation, from the Lq constraint, is the one in (4.13). Since it is 113 
that shifts by a constant we identify ex, by arguments similar to those below (2.11), 
as the boundary magnetic field and L_i as the boundary magnetization[19]. 

Finally consider W-2- If we take r(6>) = exp(9W- 2 )r then the corrections to 
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the constraints (5.2) can be computed perturbatively in 9. We find: 

d 

LX9 = L — 26>— 

LA#i = Li - + 8^A2L_ 2 L_! + • • • 

LX9 2 = L 2 - 12#A2L_iA2 + • • • 
LX9 3 = L 3 - 16#A2L_i + • • • 

= W- 2 - — 

39 

WX9 = W - 49L_ X X2 - 89X2(L_ 2 W_ 2 + 2L_ 3 W r _i) + • • • 
WX9! =Wx- 49L-! - 6#A2(L_iW_2 + 2L_ 2 W_i) + • • • 
WX9 2 = W 2 - 89X2(3W- 2 + 2L_i^_i) + • • • 

The remaining (positive n) constraints receiving no corrections to order 6>A2. 
The commutation algebra is evidently even less illuminating. The corresponding 
symmetry follows from the explicit formula for W- 2 and is highly non-local. 
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Figure Captions. 



Figure 1(a): The numerical solution to the Ising model [P, Q] = Q string 
equation (3.14) for the string susceptibility p(z), with asymptotics p(oo) = 
zXl/3 and p(— oo) = 0. The integration was performed on a mesh of 1600 
points over a range z = ±200, with a maximum error of 10A— 5. Here, z is 
the cosmological constant. 

Figure 1(b): The numerical solution to the Lee— Yang model [P, Q] =Q 
string equation (equation (2.6) with = 16/355A/C3) for p(z), with asymp- 
totics p(oo) = zXl/3 and p(— oo) = 0. An identical integration procedure to 
that used for the Ising model string equation was employed. 
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